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Abstract 

We consider first heavy particles with masses M ~ Mqjjt which arise in the Uni¬ 
verse dnring phase transitions at the Grand Unihcation scale. Using statistical mechanics 
approach we show that they behave like ideal qnantnm degenerate Bose gas which has 
the temperatnre of the Hawking thermal radiation dne to presence of the de Sitter event 
horizon. The eqnation of state for both scalar and gange bosons is presented inclnding the 
conpling constant at the Grand Unihcation scale. 


BAGS nnmbers: 04.20.Jb; 04.20.Gv; 97.60.Lf; 98.80.Gq 


According to Grand Unihed theories, hrst phase transitions occnr at the Grand Unih¬ 
cation scale Eqjjt ~ (10^^ — lfA^)GeV, and the hrst particles which are scalar and gange 
bosons get masses of order M ~ Mgut via spontaneons symmetry breaking (SSB). In 
the context of inhationary models related to a symmetry-breaking phase transition, scalar 
held reaches the minimnm of its effective potential at the end of inhation. A vacnnm en¬ 
ergy then exists in the form of coherent oscillations of a scalar held aronnd its minimnm, 
corresponding to a condensate of zero-moment um cf) particles with masses M ~ Mgut 
In models with chaotic inhation coherent oscillations of scalar held are described by a dnst 
eqnation of state corresponding to massive particles at rest 

In the present paper we consider hrst heavy nonrelativistic particles nsing statistical 
mechanics as the model-free tool. Snch an approach allows to reveal characteristic featnres 
of their behavionr which wonld not depend crncially on details of particnlar inhationary 
scenarios. 

The hrst phase transition can be related to the global discrete SSB, becanse the arrow 
of time is expected to appear then Hence scalar bosons can arise hrst in the Universe. 
Heavy nonrelativistic gange bosons are expected to emerge next via the local SSB. 
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These particles seem to be really first particles in the Universe, in accordance with 
the Tnrner ”No-Hair Theorem”: ’’Inflation greatly lessens dependence of the present state 
of the Universe npon its initial state” 

This statement goes back to the Sakharov’s ^ and diner’s ^ snggestion, formnlated 
in 1965, that the vacnnm state 

p^-e ( 1 ) 

can be the limiting state at snperhigh densities, and to the diner’s 1970 snggestion ® that 
it can be an initial state for the expanding Universe. The cosmological model describing 
evolntion of the Universe starting from the initial vacnnm state (1) has been calcnlated in 
1975. Using simple phenomenological eqnation of state for transition from an inflationary 
stage (1) to radiation-dominated stage p = e/3, it has been shown that the de Sitter 
vacnnm state (1) provides both reason for expansion and isotropy and homogeneity of the 
Universe, independently on mechanisms responsible for its arising 

In 1981 Gnth has fonnd that the eqnation of state (1) can arise in Grand Unihed 
Theories as the eqnation of state for a scalar held in a metastable state called the false 
vacnnm Since 1981 most attention was focnsed on constrncting physical mechanisms 
driving inhation and providing fnrther reheating of the Universe (for cnrrent statns of 
reheating problem see ^ and references there). 

On the other hand, in the context of SSB, the very hrst postinhationary material can 
be made np from heavy nonrelativistic particles with masses M ~ Mqut- 

Onr aim is to asses its characteristic properties which wonld be connected with main 
physical featnres of inhation and wonld not depend on details of particnlar models. 

The energy scale of hrst SSB phase transitions is U ~ Egut- At the beginning and 
dnring phase transitions the Universe is expanding in the snperlnminal way Therefore 
hrst heavy particles emerge as test particles in the vacnnm (1) with the density 

[ Egut V 

Pvac ~ PGUT — PPlanck I 

Planck/ 


They can interact within the cansally connected region conhned by the de Sitter event 
horizon 


\^TlGpyacJ \8'^/ \EgUT J 

This can be the reason why heavy particles arise hrst. Their Gompton wavelength 


( 3 ) 


A 


h 


( Epi 


Me \Egut 


t-pi 


( 4 ) 


can not exceed radins of cansally connected region (3). 

The lower limit for the particles kinetic energy follows from the nncertainty relation 


as 


P^ ^ ^ / Egut 
Mgut 3 \ Epi 


2 

Egut, 


( 5 ) 
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where p ~ fi/aQ. For Equt ~ IQ^^GeV we have 

Pgut ~ 10^^(7/cm^; ag ~ 10“^^cm; A ~ 10“^^cm; -Emin ~ lO^GeF (6) 

We can estimate the upper limit for Ekin taking into account that newborn particles will 
immediately ’’see” the Hawking thermal radiation due to presence of the de Sitter event 
horizon. Its temperature is determined by the Hawking formula 

kTff = ^ ^ ^ ^^^ ^Egut ~ 10~'^Egut ~ lO^^GeV, (7) 

where k is the Boltzmann constant. 

Since the upper limit for kinetic energy E^ax ~ kTn Egut, one can not expect 
that particles have large momenta. Therefore the CM energy of their collisions ^/s is close 
to the sum of rest energies of the colliding particles ^/s ~ 2E. As a result the dominating 
interactions are due to elastic 2 —2 body scatterings. 

For the hrst scalar particles, which are selhnteracting, we can estimate characteristic 
radius of interaction. Since particles mediating interaction are the same scalars with masses 
Msc, the cross section depending on the coupling constant Qsc is determined by 

2 

*^SC 2\2 /o\ 

C^int ~ ~ A , (SJ 


where 

The characteristic radius of interaction (dehned by the relation aint ~ ’’Lt) 

'f'int ~ CtscA ( 9 ) 

is smaller than A, if age is smaller than unity. For heavy gauge bosons (X,Y) emerging via 
the local SSB, with masses ~ MguTi the coupling constant is estimated as ^^^2 

= 9x^/471 ~ 1/40 - 1/25 (10) 

Thus, for heavy gauge bosons Tint < A due to (10). It is rather natural to expect that 
ctsc ~ Gx (in SUSY case Qsc — Qx), then the cross-section (8) is small enough for particle- 
particle interaction to produce a temperature exceeding (7). 

With the upper limit on kinetic energy given by (7), we obtain the ratio 

Ekin , E^nax f \ ( GgUt\ „_4 

This quantity characterizes the accuracy with which we can consider hrst heavy particles 
as a dust with the equation of state p = 0. 
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In order to proceed further we should make sure that we have ensemble of particles to 
apply statistics. This is just the moment when we need cosmological model of transition. 

There is no generally accepted detailed theory now, describing a transition from the 
inflationary stage of the Universe evolution to a postinflationary stage (see, e.g., ^). To 
make our estimates as model-free as possible, we describe a transition from the vacuum 
state (1) to the dust stage, using the simple phenomenological equation of state which 
is not related to a particular inflationary model and contains least possible number of 
arbitrary parameters (it has been used in 1975 to connect the inflationary onset of the 
Universe with the standard scenario of its evolution afterwards). It has the form ^^3 


p £ = £l 


(gp - 
(£0 -gl)^ 


( 12 ) 


At £ = £o it gives the vacuum equation of state (1). At £ = £i it gives the dust equation 
of state p = 0. A parameter 0 < /3 < 1 characterizes the rate of a transition. 

Integrating the Friedmann equations with the transitional equation (12), we can de¬ 
scribe evolution of a scale factor and energy density in the analytic form in the cases /3 —1, 
/9 1 and /? = 1/2. For /? = 1/2, we obtain the scale factor 


a = aoexp 



(13a) 


where A = 2a/(1 — eij 80 )! (3£i/£o). Evolution of energy density during a transition is 
given by 

£ = £0 (1 - sin^ -^) (136) 

In the case /3 -C 1 the scale factor evolves as 

a = aoexpf ^(1 -f- ^Je/eo)—] (14a) 

V2 ao/ 

and energy density as 

£ = £ 0 (^ 3 —--1^ (146) 

To estimate the duration of transition ti, we use the assumption that the whole (or 
almost whole) rest mass has been created during phase transitions at the Grand Unihcation 
scale in the form of heavy particles with masses M ~ Mgut, which later have decayed 
into some more light species, but no (or almost no) new rest mass has been created during 
reheating and during further expansion of the Universe. This assumption looks natural 
because in the course of reheating ultrarelativistic particles can be created whose total 
energy in any case can not exceed the characteristic energy E ~ Egut and whose rest 
mass energy would be therefore much smaller than this quantity. 

Applying condition of conservation of the rest mass afpi = a^odayPtoday: we obtain 

Pi ~ (1/65)pg[7t; ai ~ 10“®cm; ~ 2 x 10“^^s (15) 
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for Ptoday rsj 10 “^^g/cw?. These estimates do not depend crncially on the choice of the 
valne of the parameter /3 in the eqnation (12). The most noticeable difference takes place 
for the parameter pi which aX (3 takes the valne pi ~ {1/12Q)pqut- (The case (3^1 
corresponds to the power-law inflation.) 

The transition into the dnst implies that the whole rest mass (Mi ~ 2 X is 

contained in particles with mass Mgut ~ 2 x 10“^5r. The nnmber of particles is ~ 10^®. 
The volnme at the end of transition is Vi ~ 4 x 10“^®cm^, hence density of particles is 


n = 


Vi 




(16) 


Using the Friedmann eqnations and formnlae (13)-(14), we can check total energy bal¬ 
ance and hnd that only abont 1/4 of vacnnm energy involved in transition has been spent. 
The snbstantial amonnt remains available for farther expansion and reheating. Therefore 
the Universe is still mainly vacnnm dominated with p^ac ~ 0-7pgut, and particles can 
interact within the cansally connected regions conhned by the de Sitter horizon (3). The 
nnmber of particles in snch a region is ~ 10^^. Hence, we have an ensemble and can nse 
statistics. 

According to (16), the mean distance between particles is 

Tav ~ 10-^^cm (17) 


The characteristic radins of interaction (9) is mnch smaller than Vav and smaller than A 
(for cts < 1 and ax < 1). Therefore the free path I is determined by geometrical cross 
section ao ~ ttA^ (in hard balls approximation with A as characteristic size of particles) as 

I ~ ^ ~ 3 X 10“^^cm (18) 

As a resnlt onr particles satisfy two general criteria for ideal gas: 

< 1 Tint < I 


On the other hand, I ag which means that in cansally connected regions particles can 
collide, and their eqnation of state wonld differ from dnst eqnation of state p = 0. 

Using the npper limit for kinetic energy pnt by (7), we can estimate the lower limit 
for the de Broglie wavelength as 

H 

^deBroglie ~ ~ 10 C?7l (19) 

PH 

Since XdeBrogiie > ^av: onr ideal gas is essentially qnantnm one. Nnmber of particles is 
mnch higher than the nnmber of qnantnm states they can popnlate: 


^deBroglie A 
Tav J 


Ui 




iVi 




n- 


P^ (Uip3)/h3 {MkTnf/^ 


> 1 


( 20 ) 
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The temperature of degeneration for ideal quantum Bose gas is determined by 


14 


( 21 ) 

s Jt M 

Here Qst = 2s + 1 is statistical weight, with s = 0 for scalar and s = 1 for vector bosons. 

For our gas kTdeg ~ Because it interacts with the Hawking thermal radia¬ 

tion its temperature could be close to Th- Let us show that this is the case. Interaction 
of thermal radiation (7) with nonrelativistic heavy small particles would be essentially the 
Thomson scattering, because the condition for low-freqency limit 

Ekin ku (kTu V ^ 1 
Mc2 Mc2 "" V E J 


is satished. In our case the Thomson cross section can be written as 


StT StT 2 ^ 2 

= T W = T“ 


( 22 ) 


with a for both age and a^- The time for achieving thermal equilibrium between particles 
and radiation is the relaxation time 


T-P 7 _ 
^rel 


1 

naxc 


~ 5 X 


(23) 


Because thermal equilibrium can be achieved during a transition, and particles 

can acquire the Hawking temperature Th- Since Th T^eg, we have ideal quantum 
degenerate Bose gas with the equation of state 


P = 0.0851gst—^{kTH)^/^ (24) 

n 

We see that the process of emerging of particles looks like evaporation of a Bose condensate. 

Phenomenological cosmological model (12) we use, describes a transition which has 
hxed both beginning and the end (p = 0 at the certain moment ti). The same result has 
been obtained using smooth model with p = {■yit) — 1)£, where 7 (t) —1 at t —cx) 
Therefore this result seems to be model-independent. 

According to the standard approach (with selhnteraction A = g^^) for giving 
masses to heavy scalars via the global SSB 142 , 16 ^ have 


Mgc ~ 9scVi 

where vj is a vacuum expectation value which is expected to be vj ~ Mgut- 
The local SSB is expected to give mass to the gauge bosons of the order 


Mx ~ gxVii 
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Since both effects occur at similar energy scale E ~ Equt, one can expect that vjj ~ 
Mgut- (Note, that in SUSY case vi — vn and Qs = Qx-) Hence the equation of state for 
both scalar and gauge bosons would be the same up to the constant: 


P = const 


3/2 

{sc,x)9^sc,x)PGUT' 


2 / Egut \ 
\ Epi J 


5/2 


(25) 


The estimates presented here do not depend crucially on the choice of GUT scale as 
lO^^GeU. The scale lO^^GeU can be chosen (which is preferable by SUSY models). It is 
also possible to have both phase transitions at the same time 

We can also estimate characteristic relaxation time for particle-particle (pp) collisions 
during a transition and compare it with the duration of a transition ti. This relaxation 
time can be determined as 

pp _ ^ 

' rpJ 5 

ndintVT 

where vt is an average thermal velocity of particles. According to (7), 



Vp 



Egut 

Epi 


1/2 

~ O.Olc 


(27) 


Then 


.pp ^ ( Epi 

"" na^X^c\EGUT 


1/2 


(28) 


It represents, in particular, estimate of characteristic timescale for annihilation. We see 
that 3> ti- Therefore equilibrium state for particles cannot be achieved during a 
transition. It could mean that appropriate situation would exist at the energy scale E ~ 
Egut for generating baryon asymmetry by the mechanism proposed by Sakharov ^ which 
implies nonstationary conditions in the absence of local thermodynamic equilibrium. 

Calculating cosmological model with the equation of state (25) depending on GUT 
parameters, one will have a chance to extract some information concerning GUT physics 
from observational cosmological data, and compare it with accelerator’s data. It would 
provide a way toward astronomical testing of GUT models. 
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